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Hough Transform 



An edge is not a line... 



How can we detect lines ? 




Edges vs Boundaries 


Edges 

- Local intensity discontinuities 

- Points 

- Not dependent on models 

Boundaries 

- Composed of many points 

- May be dependent on model 



Finding lines in an image 


Option 1: 

- Search for the line at every possible 
position/orientation 

- What is the cost of this operation? 

Option 2: 

- Use a voting scheme: Hough transform 



Finding lines via Hough Transform 

• A line has a two parameters (m,b) 


Finding lines via Hough Transform 

• A line has a two parameters (m,b) 
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• Given a point (xO,yO), the lines that could pass through this 
points are all (m,b) satisfying 


Vo = m xq + b 


The equation b « >x0 m yO is a 
line in (m,b) space 


Hough Transform (Continued) 


• All points on a line in image space, yield lines in 
parameter space which intersect at a common point 

- this point is the (m,b) of the line in image space. 
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Hough Transform Algorithm 


• Initialize an accumulator array A(m,b) to zero 





Hough Transform Algorithm 


Initialize an accumulator array A(m,b) to zero 

For each edge element (x,y), increment all cells that satisfy 
b = -X m + y 
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Hough Transform Algorithm 


Initialize an accumulator array A(m,b) to zero 

For each edge element (x,y), increment all cells that satisfy 
b = -X m + y 


“This representation fails in 
case of vertical lines” 
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Polar Coordinate Representation of Line 

• p = X cos 0 + y sin 0 

- d is the perpendicular distance from the line to the origin 

- 0 is the angle this perpendicular makes with the x axis 

• Avoids infinite slope 
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Polar Coordinate Representation of Line 

• p = X cos 0 + y sin 0 


- d is the perpendicular distance from the line to the origin 

- 0 is the angle this perpendicular makes with the x axis 

• Avoids infinite slope 



The parameter space 
transform of a point is a 
sinusoidal curve 



A(p,e) 


Basic Hough transform algorithm 


1. Quantize the Hough Transform space: identify the 
maximum and minimum values of p and 0 

2. Generate an accumulator array H[d, 0]; set all values to 
zero 

3. for each edge point l[x,y] in the image 

for 0 = 0 to 180 

d = xcosO + ysinO 

H[d, 0] += 1 

4. Find the value(s) of (d, 0) \A/here H[d, 0] is maximum 

5. Calculate the equation of the line d = xcosO + ysinO 



Extensions 


Extension 1 : Use the image gradient 

• To reduce the computational load use Gradient 
information 

In step 2 : for each edge point l[x,y] in the image 

compute unique (d, 0) based on image gradient at (x,y) 

H[d, 0] +=1 

To reduce the effect of noise more than one element 
(elements in a neighborhood) in the accumulator 
array are increased 



Extensions 


Extension 2 : 

• The same procedure can be used with circles, squares, 
any other shape 



Hough Transform for Detection of 

Circles 

The parametric equation of the circle can be written as 

(x-af + iy-bf^r^ 

The equation has three parameters - a, b, r 

The curve obtained in the Hough Transform space for each 
edge point will be a right circular cone 

Point of intersection of the cones gives the parameters a, b, r 



Hough Transform for Detection of 

Circles 


Example- circles (xO,yO,r) 


(x-x0)^2 +(y-y0)^2 = r^2 


for all X 
for all y 

if E(x,y) 

for all xO 

for all yO 

r=root ( (x-xO) '^ 2 + (y-yO) '^ 2 ) 
increment H at (xO,yO,r) 

end 

end 



r 


end 


end 


end 


Detection of circle by Hough 
Transform - example 



Original Image 


Circles detected by Canny Edge 
Detector 




Detection of circle by Hough 
Transform - contd 



Hough Transform of the edge detected image Detected Circles 






